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[•"t ' Abstract 

' We analyze the time evolution of a one-dimensional quantum system with an attractive delta function potential 

whose strength is subjected to a time periodic (zero mean) parametric variation r)(t). We show that for generic 
rj(t), which includes the sum of any finite number of harmonics, the system, started in a bound state will get fully 
ionized as t — > oo irrespective of the magnitude or frequency of rj(t). For the case rj{t) — rsin(tjt) we find an 
explicit representation of the probability of ionization. There are however exceptional, very non-generic r)(t), that 
■ do not lead to full ionization. These include rather simple explicit periodic rj(t) for which the system evolves to a 

00 ' nontrivial localized stationary state related to eigenfunctions of the Floquet operator. 

m ! 

Analyse non-perturbative d'un systeme quantique modele avec force exterieure 

periodique 

\Q ' Resume 

o : 

Nous analysons revolution dans le temps d'un systeme unidimensionel avec un potentiel attractif de type 
fonction delta sounds a une variation periodique de moyenne nulle, r)(t). Nous demontrons que pour rj generique 
(en particulier pour une somme finie d'oscillations harmoniques) le systeme qui est d'abord dans un etat lie vat 
etre completement ionise pour t — > oo. Des fonctions rj(t) tres nongeneriques, toutefois explicites, existent pour 
lesquelles le systeme evolue vers un etat localise non-trivial, lie aux fonctions propres de l'operateur de Floquet 
^ associe. 

> : 

^ ! 1 Version frangaise abregee 

H | 

i Nous etudions rigoureusement le comportement pour t — > oo d'un systeme quantique unidimcnsionncl simple, avec 
potentiel attractif de type delta, soumis a une variation parametrique periodique. Dans des unites convenables, le 
Hamiltonien est de la forme 



d 2 

H(t) = H -2 v (t)6(x) = — -2S(x) -2 V (t)S(x) 

oil Hq a un seul etat lie it& = e - ' 21 ' d'energie uio = —1 et un spectre continu sur l'axe reel positif, avec fonctions propres 
generalisees, voir eq. ©. 

On peut developper la solution de l'equation de Schrodinger ip(x,t) par rapport aux fonctions propres de Hq ||S} 
avec conditions initielles 9(0) — do, O(fc,0) = Oo(fc) normalisees convenablement, eq. ©. Alors, la probabilite de 
survie de l'etat lie est P(t) — \6(t)\ 2 , alors que \Q(k,t)\ 2 dk donne la "fraction de particules ejectees" avec (quasi-) 
impulsion dans l'intervalle dk. En prennant la fonction Y donnee par Q on obtient les equations eq. ijHJ) et Y 
satisfait une equation integrate, QlOf) . Notre methode d'analyse utilise les proprietes analytiques de la transformation 
de Laplace y(p) de Y(t) pour determiner les proprietes asymptotiques de Y(t) par rapport a t. 

*Department of Mathematics, Hill Center, Rutgers University, New Brunswick, NJ 08903, USA. e-mail: costin@math.rutgers.edu, 
rcostin@math.rutgers.edu, lebowitz@math.rutgers.edu, rokhlenk@math.rutgers.edu. 



1 



2 Cas ou 77 (t) est harmonique 



Theoreme 1 Si rj(t) = rshiut , la probabilite de survie \9{i)\ 2 tend vers zero quand t — > 00, pour tous oj > et 

Remarques. (1) On obtient une formule exacte pour 9(t) oil F w est periodique de periode 2ttui~ 1 et les fonctions 
h m satisfont 1(12(1, Pas trop pres des resonances, si |w — n _1 | > 0(r 2_<5 ) pour tout entier positif n, |i*L(i)| = 1± 0(r 2 ) 
et les coefficients de Fourier de F decroissent plus vite que r' 2 ™ 1 ' |m| - '™ 1 '/ 2 . Aussi, la somme en (|llfl est plus petite que 
0{r 2 t~ 3 / 2 ) pour t grand, et les h m decroissent avec m plus vite que r'™ 1 '. 

(2) On voit grace a ((11(1 que pour des temps d'ordre 1/T oil V = 2di(j), la probabilite de survie pour u> pas trop 
pres d'une resonance decroit comme exp(— Tt), et que le comportement asymptotique en t est |#(i)| 2 = 0(t~ 3 ) avec 
beaucoup d'oscillations. 

(3) Quand r est plus grand, le comportement polynomial-oscillatoire commence plus tot et la probabilite de survie 
est plus grande. Ce phenomene est parfois appele stabilisation atomique. 

(4) En utilisant une fraction continue convergente, on peut calculer T pour ui and r arbitraires. Pour r petit, si n 
est la partie entiere de uj^ 1 + 1 ct si oj^ 1 ^ N, alors, pour T > (t = r^ 2n T) 1 T est donne par 

(5) Le comportement de V est different aux resonances uj^ 1 € N. Par exemple si uj — 1 = r 2 /\/2 on trouve la 
formule (JT3J. 

2.1 Cas periodique general 

On ecrit 77 sous la forme ljT5)l et nos hypotheses sur les Cj sont (a) ^ rj S L°°(T), (b) Co = et (c) C-j — Cj. 
Considerons aussi l'hypothese de genericite (g) suivante: on definit la translation a droite T(Ci, C%, C n , ...) = 
(C2, C3, C n +i, ...)■ Alors, C £ ^(N) est generique par rapport a T si l'espace de Hilbert engendre par toutes les 
translations de C contient le vecteur e\ = (1,0,0...,), cf. 1)16(1 . Un cas important est donne par les polynomes 
trigonometriques. Un exemple qui ne satisfait pas a l'hypothese (g) est ifTTll pour A € (0, 1), pour lequel C„ = — rX n 
for n > 1. 

2.2 Resultats dans le cas periodique 

Theoreme 2 [10J. Sous les hypotheses (a), (b), (c) et (g), la probabilite de survie Pit) de I'etat lie u^, |#(i)| 2 tend 
vers zero quand t — > 00 . 

Theoreme 3 [10 . Pour ipo(x) — Ub(x) il existe des valeurs de \, 10 et r en jJ7| ), pour lesquelles \0(t)\ -f^ quand 
t — > 00. 

Remarque. Le Theoreme 3 peut etre etendu pour demontrer que pour r et w donncs cn 1(17(1 il existe un ensemble 
iniini de A, avec un point d' accumulation en 1, pour lequel 8(t) -f^ 0. 

3 English version 
3.1 The problem 

We are interested in the nature of the solutions of the Schrodinger equation 

iW t ^=[H Q + H 1 {t)]i> (1) 

Here ip is the wavefunction of the system, belonging to some Hilbert space 7i, Hq and Hi are Hermitian operators and 
equation 0) is to be solved subject to some initial condition ipo. Ho has both a discrete and a continuous spectrum 
corresponding respectively to spatially localized (bound) and scattering (free) states in E rf . Starting at time zero with 
the system in a bound state and then "switching on" at t = an external potential Hi(t), we want to know the 
"probability of survival", P(t), of the bound states, at times t > 0; 1 — P(t) is the probability of ionization P3-[5|- 

When uo is sufficiently large for "one photon" ionization to take place, i.e., when hio > —Eq, Eq the energy of the 
bound (e.g. ground) state of H and r is "small enough" for Hi to be treated as a perturbation of H then the long 
time behavior of P(t) is asserted in the physics literature to be given by P(t) ~ exp[— Tpt]. The rate constant Tp is 
computed from first order perturbation theory according to Fermi's golden rule. It is proportional to the square of 



2 



the matrix element between the bound and free states, multiplied by the appropriate density of continuum states in 
the vicinity of the final state which will have energy hui + Eq [HI EI IB] ■ 

The results described here show that the phenomenon of ionization by periodic fields is very complex indeed once 
one goes beyond the perturbative regime. 



4 Our model 

We consider a very simple quantum system where we can analyze rigorously many of the phenomena expected to 
occur in more realistic systems described by JQ. This is a one dimensional system with an attractive delta function 
potential. The unperturbed Hamiltonian Hq has, in suitable units, the form 

d 2 

H = -—— - 25(x), -oo < x < oo. (2) 
ax z 

Hq has a single bound state Ub(x) — e^' 1 ' with energy — ojq = — 1. It also has continuous uniform spectrum on the 
positive real line, with generalized eigenfunctions 



u(k, x) = -= e lkx rrre* 1 ** 1 , -oo < k < oo (3) 

y/2lv V l + i\k\ J 

and energies k 2 . 

Beginning at t — 0, we apply a parametric perturbing potential, i.e. for t > we have 

H(t)=H -2r 1 (t)6(x) (4) 

and solve the time dependent Schrodinger equation Q for ip{x,t), with ip(x, 0) = ipo(x). Expanding ip in eigenstates 
of Hq we write 

if>(x,t) = 0{t)u b (x)e lt + / Q(k,t)u(k,x)e' lk l dk (t > 0) (5) 



with initial values 6(0) = 6q, Q(k,0) = &o(k) suitably normalized, 

/•oo 

(^o,^o) = |^o| 2 + / \& (k)\ 2 dk = l (6) 



We then have that the survival probability of the bound state is P(t) = ^(t)! 2 , while \Q(k, t)\ 2 dk gives the "fraction 
of ejected particles" with (quasi-) momentum in the interval dk. 

This problem can be reduced to the solution of an integral equation in a single variable |5|-|12|. Setting 

Y(t) =^(x = 0,t) V (t)e u (7) 

we have 



6(t) =6 + 2i f Q Y(s)ds , 
G(k, t) = Q (k) + 2|fc|/[V2iF(l - i\k\)] f* Y(s)e^ 1+k2 >ds 

Y(t) satisfies the integral equation 



Y 



(t) = r)(t) + J [2i + M(t ~ t')]Y(t')dt'^ = rj(t)(l(t) + (2i + M) * 



(8) 
(9) 



(10) 



where the inhomogeneous term is 



I(t) = 6 a + 



j f°° eo(k) + e Q (-k) c _ l(k 2 +1)t 

" ^ Jo 1 + ik 



dk, 
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and 



with 

f*9= / f{s)g{t- s)ds 
Jo 

Outline of the technical strategy. The method of analysis in |5J^], relies on the properties of the Laplace 
transform of Y, y(p) = CY{p) = e- pt Y(t)dt. 

We show that I jlOjl has a unique solution in suitable norms. This solution is Laplace transformable and the Laplace 
transform y satisfies a linear functional equation. The solution of the functional equation is unique in the right half 
plane provided it satisfies the additional property that y{po + is) is square integrable in s for any po > 0- We use the 
functional equation to determine the analytic properties of y(p). 

This is done using (appropriately refined versions of) the Fredholm alternative. After some transformations, 
the functional equation reduces to a linear inhomogeneous recurrence equation in 1%, involving a compact operator 
depending parametrically on p. The dependence is analytic except for a finite set of poles and square-root branch- 
points on the imaginary axis and we show that the associated homogeneous equation has no nontrivial solution. We 
then show that the poles in the coefficients do not create poles of y, while the branch points are inherited by y. The 
decay of y(p) when |9(p)| — > oo, and the degree of regularity on the imaginary axis give us the needed information 
about the decay of Y(t) for large t. 

4.1 Case when rj(t) is harmonic 

Theorem 1 (|_9|) When 77(f) — rsinwf the survival probability \0(t)\ 2 tends to zero as t — ► oo, for all lo > and 

Remarks. (1) The detailed behavior of the system as a function of t, cu, and r is obtained from the singularities of 
y(p) in the complex p-plane. We summarize them for small r; below 6 > 0. For definiteness we assume in the following 
that r > 0. 

At p — {inuj — i : n £ Z}, y has square root branch points and y is analytic in the right half plane and also 
in an open neighborhood M of the imaginary axis with cuts through the branch points. As — > oo in Af we 
have \y(p)\ = 0(ru\p\~ 2 ). If \ui — i| > const n O(r 2 ~' 5 ), n e Z + , then for small r the function y has a unique pole 
Pm — Po + imuj in each of the strips — muj > ^s(p) + 1 ± 0{r 2 ~ s ) > —mu — cj, m € Z. dt(p m ) is strictly independent 
of to and gives the exponential decay of 8. Laplace transform techniques show that 



6(t)=e-^'> t Fjt)+ e. (miu> - i]t h m {t) (11) 

in— — oo 

where F u is periodic of period 2ttuj~~ 1 and 

oo 

M*) ~ £ c mJ t~ 3/2 ~ j as t -» oo, ar S(*) e ( - I - 0, 1 + 0) (12) 

Not too close to resonances, i.e. when \uj — 77 _1 | > 0(r 2 ~ s ), S > 0, for all integer n, \F u (t)\ = l±0(r 2 ) and its Fourier 
coefficients decay faster than r' 21 ™' |m|~' m '/ 2 . Also, the sum in does not exceed 0(r 2 t -3 / 2 ) for large t, and the h m 
decrease with m faster than r' m '. 

(2) By jn) , for times of order 1/r where V = 23?(7), the survival probability for u> not close to a resonance decays 
as exp(— Tt). It follows from our analysis that for small r the final asymptotic behavior for t — > oo, is |6>(i)| 2 = 0(t~ 3 ) 
with many oscillations as described by 

(3) When r is larger the polynomial-oscillatory behavior starts sooner. Since for small r the amplitude of the 
late asymptotic terms is 0(r 2 ), increased r yields higher late time survival probability. This phenomenon, sometimes 
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referred to as atomic stabilization, can be associated with the perturbation-induced probability of back-transitions to 
the well. 

(4) Using a continued fraction representation T can be calculated convergently for any to and r. The limiting 
behavior for small r of the exponent T is described as follows. Let n be the integer part of uj^ 1 + 1 and assume 
lu^ 1 £ N. Then we have, for T > (t = r~ 2n T), 



f = -T _1 lim In \6(r~ 2n T)\ 2 = v = (13) 



)-2n+2 



V nu) — 1 



i — 



nu 



Y[ (i - vi - ™) ! 



(5) The behavior of T is different at the resonances to 1 £ N. For instance, whereas if ui is not close to 1, the 
scaling of T implied by 113(1 is r 2 when w > 1 and r 4 when | < w < 1, by taking ui — 1 = r 2 /V2 we find 

2 i/4 3 / 4 

-T- 1 lim In 6»(r~ 3 r) = (14) 

r^o 1 v ;| 8 16 v ; 

w=l+r 2 /\/2 

4.2 General periodic case 

We write 



J2 (C j e iuit + C-je - ^*) (15) 



»7 = 

i=o 

Our assumptions on the Cj are (a) ^ ?] £ L°°(T), (b) Co = and (c) C_j = Cj. 
Genericity condition (g). 

Consider the right shift operator T on ^(N) given by T(Ci, C2, C n , ...) = (C2, C3, C n +i, ...). We say that 
C 6 ?2(N) is generic with respect to T if the Hilbert space generated by all the translates of C contains the vector 
ei = (1,0,0...,) (which is the kernel of T): 

OO 

e a e Y T n C (16) 

n=0 

(where the right side of ((16(1 denotes the closure of the space generated by the T n C with n > 0.) This condition is 
generically satisfied, and is obviously weaker than the "cyclicity" condition h(?$) G \J^ = oT n C = {0}, which is also 
generic. 

An important case, which satisfies 1(161) . (but fails the cyclicity condition) corresponds to 77 being a trigonometric 
polynomial, namely C ^ but C n = for all large enough n. 
A simple example which fails 1 (16(1 is 

„ , A — cos(ojt) 
/w l + A 2 -2Acos(wi) v ; 

for some A € (0, 1), in which case C n — —r\ n for n > 1. In this case the space generated by T™C is one-dimensional. 
We prove that there are values of r and A for which the ionization is incomplete, i.e. 9 it) does not go to zero for large 
t. 

4.3 Results in the periodic case 

Theorem 2 ([10J) Under assumptions (a), (b), (c) and (g), the survival probability Pit) of the bound state it& , |0(i)| 2 
tends to zero as t — > 00. 

Theorem 3 ([10j) For ipo{x) — Ub{x) there exist values of \, uj and r in jlTj) , for which \6(t)\ -/-* as t — ► 00. 
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Remarks. 

1. Theorem [2 can be extended to show that J D \tp(x,t)\ 2 dx — > for any compact interval D S R. This means that 
the initially localized particle wanders off to infinity since by unitarity of the evolution J R \tp(x,t)\ 2 dx — 1. Theorem^ 
can be extended to show that for some fixed r and ui in l|17|l there are infinitely many A, accumulating at 1, for which 
0(f) -f* 0. In these cases, it can also be shown that for large t, 9 approaches a quasiperiodic function. 

2. There is a direct connection between our results and Floquet theory where, for a time-periodic Hamiltonian 
H(t) with period T — 2n/u>, one constructs a quasienergy operator (QEO). Complete ionization thus corresponds to 
the absence of a discrete spectrum of the QEO and conversely stabilization implies the existence of such a discrete 
spectrum. In fact, an extension of Theorem [3] shows that for the initial condition ipg — u&, tpt approaches such a 
function with /j, = —so. More details about Floquet theory and stability can be found in |14M15| 

3. We are currently investigating extensions of our results to the case where Hq = — V 2 + Vq(x), x S M d , has a 
finite number of bound states and the perturbation is of the form 77(f) V\ (x) and both Vq and V\ have compact support. 
Preliminary results indicate that, with much labor, we shall be able to generalize Theorem|21 to generic V\{x). The 
definition of genericity will, however, depend strongly on Vq. 
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